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Summary. After a short introduction to the characteristic geometry underlying weakly 
hyperbolic systems of partial differential equations we review the notion of symmetric 
hyperbolicity of first-order systems and that of regular hyperbolicity of second-order sys- 
tems. Numerous examples are provided, mainly taken from nonrelativistic and relativistic 
continuum mechanics. 
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' 1 Introduction 

o 

The notion of hyperbolicity of a partial differential equation (PDE) , or a system of 
■ . PDE's, is central for the field theories of mathematical physics. It is closely related 

I to the well-posedness of the Cauchy problem and to the causal structure underlying 

these theories. In standard theories describing relativistic fields in vacuo this causal 
structure is that given by the spacetime metric, a second-order symmetric tensor of 
Lorentzian signature. If matter is included, things become both more complicated 
and more subtle. In fact, the awareness of some of those complications predates 
Relativity by centuries. An example is afforded by the phenomenon, already studied 
by Huygens, of birefringence in crystal optics ^ . 

There is currently an increase of attention in the field of Relativity, due in part 
^ ' to demands from Numerical Relativity, devoted to certain notions of hyperbolicity 

applied to the Einstein equations (for an excellent review see [18]). There the main 
challenge, not discussed in the present notes at all, comes from the fact that, already 
in vacuum, the Einstein equations by themselves, i.e. prior to the imposition of any 
gauge conditions, are not hyperbolic. The main burden, then, is to find a "hyperbolic 
reduction" turning the Einstein equations, or a subset thereof, into a hyperbolic 
system appropriate for the purpose at hand. However the complications in the 
causal structure one finds in continuum mechanics, which are our main focus here, 
are absent in the Einstein vacuum case - at least for the reductions proposed so far. 
Of course, these complications do come into play ultimately once matter-couplings 
are included. 

These notes attempt an elementary introduction to some notions of hyperbolicity 
and the "characteristic geometry" associated with or underlying these notions. The 
section following this one is devoted to the general notion of a hyperbolic polyno- 
mial, which in our case of course arises as the characteristic polynomial of a PDE. 
It is interesting that this notion is on one hand restrictive enough to encode essen- 
tially all the required features of a theory in order to be "causal" - on the other 
hand flexible enough to account for an amazing variety of phenomena - relativistic 
or nonrelativistic - ranging from gravitational radiation to water waves or phonons 

^ For a fascinating account of the history of the associated mathematics see [20]. 
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in a crystal. We devote a significant fraction of Sect. 2 to examples, which at least in 
their nonrclativistic guise all appear in the standard literature such as [13], though 
not perhaps from the unified viewpoint pursued here. Some of these examples are 
not fully worked out, but perhaps the interested reader is encouraged to fill in more 
details, possibly using some of the cited literature. We hope that some workers in 
Relativity, even if they have little interest in continuum mechanics for its own sake, 
find these examples useful for their understanding of the notion of hyperbolicity. 
While hyperbolicity of the characteristic polynomial of a theory is important, it 
is not in general sufficient for the well-posedness of the initial value problem for 
that theory. Well-posedness is the subject of our Sect. 3. We recall the notion of a 
symmetric hyperbolic system of a system of 1st order PDE's, which is indeed suffi- 
cient for well-posedness. A similar role for 2nd order equations is played by a class 
of systems, which were to some extent implicit in the literature, and for which an 
elaborate theory has been recently developed in [10, 11]. These systems are called 
regular hyperbolic. They encompass many second order systems arising in physics 
one would like to qualify as being hyperbolic - such as the Einstein equations in 
the harmonic gauge. If applicable, the notion of regular hyperbolicity is particularly 
natural for systems of 2nd order derivable from an action principle, as is the case 
for many problems of continuum mechanics. We show the fact, obvious for symmet- 
ric hyperbolic systems and easy-to-sce although not completely trivial for regular 
hyperbolic ones, that these systems are special cases of weakly hyperbolic systems, 
i.e. ones the determinant of whose principal symbol is a hyperbolic polynomial. We 
also touch the question of whether a system of the latter type can be reduced to one 
of the former type by increasing the number of dependent variables. Throughout 
this section our treatment will be informal in the sense of ignoring specific differen- 
tiability requirements. We also do not touch questions of global wellposedness. 

2 Hyperbolic Polynomials 

The PDE's we are interested in are of the form 

M^--'" {x, f,df,..., . . . + lower order terms. (1) 

Here A,B = 1, ...m and fii = 1, . . .n. Relevant equations of this form are the Euler 
equations for a barotropic fluid (for n = 4, ^ = 1, m = 4), the Einstein equations 
(for n — A^l = 2,m = 10) or the equations governing an ideal elastic solid (for 
n = 4,1 = 2,m = 3). The Maxwell equations, in the form they are originally 
written down, are not of this form, but a suitable subset of them is, as we will 
discuss later. 

The principal symbol of the PDE (1) is defined as 

MAB{k) = M^^s-^^k^,...k^„ k^G{R"r (2) 

We here suppress the dependence on x and on /. The characteristic polynomial P{k) 
is defined by P{k) = det MAB{k), where the determinant is taken with respect to 
some volume form on /-space C M™. The polynomial P{k) is homogenous of degree 
p = m-l. A homogenous polynomial of degree p > is called hyperbolic with respect 
to e (M")* if P(^) 7^ and the map A ^ P{r] + A^), itself a polynomial of degree 
p, has only real roots Aj, i = 1, • ■ -p for all rj € (R")*. 

The roots Xi{^, rj) need not be distinct. If, for all ry with rjA^ ^0, Xi{^, rf) ^ Xj (^, t]) 
for i ^ j, P is called strictly hyperbolic ^). We write G* for the set of fc G (R")* \{0}, 

^ This case is not general enough for the purposes of physics. Furthermore there exist 
physically relevant cases of non-strictly hyperbolic polynomials which are stable, in the 
sense that they possess open neighbourhoods in the set of hyperbolic polynomials just 
containing non-strictly hyperbolic ones [26,28]. 
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where P vanishes. It is sometimes called the cone of characteristic conormals. 

It is clear that a product of hyperbolic polynomials is hyperbolic. Also, if a hyper- 
bolic polynomial can be factorized into polynomials of lower degree (in which case 
it is called reducible), these factors are also hyperbolic. There is a wealth of infor- 
mation which can be inferred about a polynomial P{k) if it is hyperbolic. Before 
explaining some of this, we look at a few examples for hyperbolic polynomials. 

Example 1: P{k) = {a,k) = a'^k^ for some nonzero a** G R". The set C* is a 
punctured hyperplane C (R")*. 

Clearly P{k) is hyperbolic with respect to any such that a^^fi ^ 0. The polyno- 
mial P{k) = (ai, fc)(a2, fc)(a3, fc), with ai,a2,a^ linearly independent G M'^, arises 
in the problem of finding, for a three dimensional positive definite metric, a coordi- 
nate system in which the metric is diagonal (see [14]) - which shows that hyperbolic 
problems can also arise in purely Riemannian contexts. 

Example 2: P{k) = ^'^"k^k^, where 7'*'' is a (contravariant) metric of Lorentzian 

signature (— , +,..., +). The set 6* is the two-sheeted Minkowski light cone. 
When n = 2, P{k) is hyperbolic with respect to any non-null when n > 2, P{k) 
is hyperbolic with respect to any with "f^'^^^^v < 0, i.e. ^ is timelike with respect 
to 7^"^. Checking that P{k) is hyperbolic according to our definition is equivalent 
to the so-called reverse Cauchy-Schwarz inequality for two covectors one of which 
is timelike or null with respect to 7'"' (which is the mathematical rationale behind 
the the twin "paradox" of Relativity). Surprisingly there are similar inequalities 
for general hyperbolic polynomials (see [19]) which play a role in diverse fields of 
mathematics [5]. 

Example 2 is of course the most familiar one. If it arises from nonrelativistic field 
theory, the quantity 7^"^ currently runs under the name of the "Unruh or acoustic 
metric" [4] (see also [12]) in the Relativity community. It is not an elementary object 
of the theory, but is built as follows: Take first the Galilean metric /i'"', a symmetric 
tensor with signature (0, -|-, . . . , -|-) together with a nonzero covector satisfying 
h^^Ti, = 0: these are the absolute elements . Then pick a 4- vector normalized so 
that u^T^ = 1 and define 7^*^ = h^'^ — cr'^v/^u'^ . This describes waves propagating 
isotropically at phase velocity c in the rest system, defined by m^, of a material 
medium. The relativistic version of the above is as follows: Start with the spacetime 
metric g^'^ and define 7^*^ = g^^ + (1 ^ cr'^)u'^v^ , where u/* is normalized by taking 
Tjx — ~9iJ.vu'^ J with Qfj^i, the covariant spacetime metric defined by g^^g" = • 
Note: if there are metrics 7^ ,72"^ with C2 < Ci, then the "faster" cone lies inside 
the slower one. We will come back to this point later. 

Example 3: P{k) = s^'^k^k^,, where s^'^ has signature (— , -|-, . . . , -|-, 0, . . . , 0), is 

hyperbolic with respect to any ^ such that s^^£,^£^i, < 0. 

Here is a case occurring in the real world. Let g^;^ be a Lorentz metric on M^, 
u'^ a normalized timelike vector, i.e. g^vU^u" = —1, F^i, = F^^^] nonzero with 
FfiuW = 0. The quadratic form s''" = -eut^u" + 1/2 Fp^PP'^gi"' - F^'pF^'P, with 
e > 0, has signature (— , -|-, -|-, 0). The characteristic cone C* of P{k) = consists 
of two hyperplanes punctured at the origin. When e is interpreted as e = "energy 
density -I- pressure" and as the frozen-in magnetic field of an ideally conducting 
plasma, then P{k) describes the Alfven modes of relativistic magnetohydrodynamics 
[45] [2]. 

Example 4: Let n = 4, e^^xp some volume form on M'' and m^^^P = m'^'^ll'^^l. With 

Q^upa ^ g^^^^ g^^^^ j^al}K{^^u\5<t>\p^<T)e^u ^^gfl^g p pQ^-^ ^ G''"'"^ k^k^kpk^. 

As a special case take m'"'^P of the form mt"'^P = h^^t^h^^P - e^^Pii^^uP + cP^Pu^^u^, 
where the symmetric tensors hP" and sP" , both of signature (0, -|- -|- -|-), satisfy 
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h^'^Ti, = e^'^Ti, = for u^Tfj, ^ 0: this is the situation encountered in crystal optics 
with the nonzero eigenvalues of e^"^ relative to ft,'"^ being essentially the dielectric 
constants. The crystal is optically biaxial or triaxial, depending on the number of 
mutually different eigenvalues. The 4th order polynomial P{k) turns out to be hy- 
perbolic with respect to all in some neighbourhood of = , and the associated 
characteristic cone is the Fresnel surface (see e.g. [25]). For an optically isotropic 
medium or in vacuo P{k) is reducible, in fact the square of a quadratic polynomial 
of the type of Example 2. We leave the details as an exercise. More general condi- 
tions on m^^^f in order for P{k) to be hyperbolic can be inferred from [32]. 
The quartic polynomial P{k), as defined above, comes from a generalized ("pre- 
metric") version of electrodynamics (sec [24]), as follows: Let be the electromag- 
netic field strength and write H^^" = m^^^Txp for the electromagnetic excitation. 
The premetric Maxwell equations then take the form 

d[t,{e,x],„HP^) = J^.A, dy^F^x] = 0, (3) 

where J/^i^x is the charge three form The Eq.'s (3) reduce to the standard ones 
in vacuo when m^"^'' ~ g^\-t^g"\P with g'^" the metric of spacetime. If one sets 
^Mi'Ap ^ ^A[M^-^]p^ ^ith = h^"" - c-'^u^'u'', h^"' the Galilean metric and u^' a 
constant vector field s.th. u^^t^ = 1, one has the Maxwell equations in a "Galilean" 
(not Galilean-invariant) version with u** describing the rest system of the aether 
(see [43]). One then looks at hypersurfaces along which singularities can propagate. 
The result is that the conormal of such surfaces has to satisfy P{n) = 0. Put 
diff'erently, one can look at the 8x6- principal symbol of the Maxwell equations: 
then P(fc) = is exactly the condition for this principal symbol to have nontrivial 
kernel. If one considered an appropriately chosen subset amongst (3), the evolu- 
tion equations, one would obtain an equation of the form (1), whose characteristic 
polynomial contains P{k) as a factor. We will treat the vacuum case of this later. 

Our last and most complicated example comes from elasticity [6], namely 
Example 5: Take n = 4, Z = 2, m = 3 in Eq.(l) with 

MZ, = -GABU^'vr + C'X"B, (4) 

where Gab = G(ab) ^md = C^'^ and C^gT^ = for some covector r satisfying 
(u, r) ~ u'^Tfj_ — 1. The theory is intrinsically quasilinear: all quantities entering 
Eq.(4) are functions of / and df and in general also of x. For example /"^ is 
required to have maximal rank, and u^^ satisfies m^(9^/'^) = 0. Furthermore G'^^ = 
GADBE{^pf^){^„f^)hP^'h'"' , with h^"',T^, being, in the nonrclativistic case, the 
absolute Galilean objects, or, in the relativistic case, /i^" = g'^^ -|- m^m" and = 

There are the following basic constitutive assumptions 

Gab is positive definite G'^gin^m^ rjnTj^ > for m 7^ 0, A r 7^ (5) 

Defining the linear map (M)^b by {M)^B{k) = -{u, kf S^B + G^^G'^j^gk^.k^, the 
polynomial P{k) can, by general linear algebra, be written as 

6 P{k) = {trMf - 3 {trM^){trM) + 2 trM^. (6) 

It will follow from a more general result, to be shown below, that the 6th order 
polynomial P{k) is hyperbolic with respect to in some neighbourhood of r^. In 
the special case of an isotropic solid, the "elasticity tensor" Gabde has to be of the 
form 

^ These equations play a certain role in current searches for violations of Lorentz invari- 
ance in electrodynamics [30] 
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Cabde — I GabGde + 2m Gjj(^aGb)e 



(7) 



and the second of Eq.(5) is satisfied iff C2 = m > 0, = l + 2m > 0. The polynomial 
P{k) turns out to reduce to the form 



of pressure and shear waves respectively. If the medium is elastically anisotropic, 
such as a crystal, one can start by classifying possible fourth-rank tensors Cabde 
according to the symmetry group of the crystal lattice, allow for dislocations, etc. 
The richness of possible structure of C* and the corresponding range of captured 
physical phenomena - studied by theoreticians and experimentalists - is enormous. 

This c'nds our list of examples. We now turn to some general properties of hyper- 
bolic polynomials and their physical interpretation. It is clear from the definition 
that e* has codimension 1: since P{r] + X^,) has to have at least one complex root 
for each 77, and the roots are all real, there is at least one real root. And since there 
are no more than p ■ I different roots, C* can not have larger codimension. It is 
then known from real algebraic geometry that C* consists of smooth hypersurfaces 
outside a set of at least codimension 2 (see [9]). The roots Aj(^,T/) can for fixed ^ 
be assumed to be ordered according to Ai < A2 < . . . < Aj, for all 7]. The set of 
points k = 1] + Xi{^,r])^ is called the i'th sheet of C*. The hypersurface C* has to be 
smooth at all points k lying on a line intersecting p different sheets ^. In particular 
all sheets are everywhere smooth when P is strictly hyperbolic. 
Next recall that the defining property of a hyperbolic polynomial refers to a par- 
ticular covector ^. That covector however is not unique. It is contained in a unique 
connected, open, convex, positive cone r*{^) of covectors ^' sharing with ^ the 
property that P{^') ^ and P{ri + A^') has only real zeros Ai(^', 77) [19]. Note that 
r*(^) = -r*(-$). Furthermore dr*((^) C 6*, and r*(^) is that connected com- 
ponent of the complement of C* containing ^. Not all points of dr*{S,) have to be 
smooth points of 6*. 

The roots Xi{^,r]),i = 1, . . .p, due to the homogeneity of P, are homogenous in ^ 
of order —1 and positively homogenous of order 1 as a function of rj. They also 
satisfy Ai(^,r/) = — Ap+i_i(^, 77). At regular points of 6*, i.e. when the gradient of 
P at 77 + Aj^ is non-zero. A, 77) is a smooth function of its arguments due to the 
implicit function theorem. 

Next choose a vector X e M" so that {X,C) > for all e r*{0- We call such 
a vector "causal". Wc now look at the intersection S of the hyperplane {X,^') = 1 
with e*. Note {X, = 1 is transversal to 6* at smooth points of G*,so S is smooth 
there also. Note also that S may be empty, as in Ex.1. For reasons explained below, 
S is often called "slowness siirfacc" . To describe S more concretely, we pick some 
T e r*{^) with {X,t) = 1. The pair X,t constitutes a "rest frame". Using it, we 
can decompose every covector k as fc,^ = + A:^ where fc-"- is tangential to the 
hyperplane, i.e. (X, fc-*-) = 0. This k lies on C* iff Ai(r, fc^) = 1 for some i. Thus 
S consists of the sheets Aj(r, fc^) = 1, viewed as {n — 2)-surfaces in k^ € M""-^. 
Here J runs through some subset of the i's parametrizing A^ from above. Clearly, as 
J increases, these sheets form a nested family of not necessarily compact surfaces. 
^ The innermost of these surfaces is nothing but the intersection of dr*{£^) with 
the hyperplane {X,^') = 1 and is hence convex. We call C*(^) those components 

* The reason is that a polynomial of order p in one real variable, if it has p different zeros, 
has non-vanishing derivative at each zero, so fc is a non-critical point of P. 

^ In particular, when sheets seem to pass through each other, the two sides are counted 
as belonging to different sheets. 
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of C* which consist of half-rays connecting the origin with the points of S. In the 
examples 2, 4, 5 the set S and C*(^) consist of at most 1, respectively 2 and 3 sheets. 
For the last-mentioned case, see [15]. Not all the occurring sheets are compact. It is 
possible for example for P{k) to be an irreducible hyperbolic polynomial with some 
sheets of S compact and others non-compact: this is the case e.g. for the acoustic 
modes in magnetohydrodynamics [13]. 

We now explain the name "slowness surface" . Consider the hyperplane in M" given 
by (a;, fc) = for fixed k G G* , i.e. the wave front of the plane wave associated with k. 
To measure the "speed" at which this wave front moves, decompose observers with 
tangent V which "move with this wave front" , i.e. such that (V, k) = 0, according 
to V = X + v-^. It follows that (v-*-, fc-*-) = —1. Thus, if there is a natural "spatial" 
metric h mapping elements l-^ into elements w-^ = hol-^ orthogonal to r, one can 
define the "phase velocity" wj^j^ = — h o fc-'-. Thus, the smaller fc-*- is, the 
larger the phase velocity. Of course the equation (t;-*-, fc-"-) = — 1 does not define v-^ 
uniquely. But there is a "correct" choice for V tangential to the wave front, called 
"ray or group velocity", which is independent of any spatial metric, and which is 
defined at least when fc is a smooth point of C*: this V is given by the conormal 
to e* C (K")* at fc, which by duality is a vector e K". If k is in addition a non- 
critical point, this ray velocity is ~ 9/9fc^P(fc), which satisfies fc^y = by 
the positive homogeneity of P{k). The spatial group velocity in the frame X, r can 
then be written as 



which is also the textbook expression. 

We should add a cautionary remark here. Although the differential topology of the 
slowness surface is independent of the choice of X satisfying {X,^') > for all 
^/ g r*{^), its detailed appearance, and physical quantities such as phase velocity, 
group velocity or angle between two rays do of course depend on the choice of 
rest system X, t and a notion of spatial metric with respect to that observer. Of 
course there will be, for any particular physical theory, a singled-out class of rest 
systems, e.g. can be the absolute object in a Galilean spacctime or be of the 
form = —gij.vX" in a relativistic theory. Or the slowness surface can have more 
symmetry (say symmetry with respect to reflection at the origin) in some rest 
system than in others, as is the case with crystal optics or elasticity. For a careful 
discussion of these issues, in the more specialized context of "ray-optical structures" 
on a Lorentzian spacetime, consult [34] . 

We now come back to the "ray" concept. If fc is a smooth critical point of 6*, finding 
the map k V{k) is already a nontrivial problem in algebraic geometry [35]. If fc 
is not a smooth point of C*, there is no unique assignment of a group velocity to k. 
Still well-defined is the set C of all ^ 7^ satisfying 



called the dual or ray cone. Loosely speaking, each sheet of the ray cone corresponds 

to a spherical wave front tangent to (or "supported by") the planar wave fronts 
defined by the different points k in some corresponding sheet of C* [13]. There holds 
(e*)* = e. The dual cone is again an algebraic cone, which, except in degenerate 
cases, is again the zero-set of a single homogenous polynomial. The structure of this 
dual cone, in particular its singularity structure which can be very complicated, 
is another difficult matter of real algebraic geometry. For example the degree of 
its defining polynomial is in general much higher than that of C* (see [37], [21]). 
This "dual" polynomial need not be hyperbolic: in order to be hyperbolic it would 





(V, fc) = where P(fc) = 0, 



(10) 
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have to have a central sheet which is convex, which is not the case for some of the 
examples one finds in the literature. In our examples from above the situation is 
as follows: In our Ex.1 the dual cone C* consist of the two half-lines {aa'^|a > 0} 
and {aa^\a < 0}. The cone dual to the quadratic cone g'^^kfjikv = in Ex.2 is 
given by g^uV'^V — with g^i^g"^ = S^"^. For a nonrelativistic acoustic cone 
'y^'^ = h^" — {l/c^)u^u'' we obtain for the ray cone "^^v = h^v — c^t^Tu, where 
hiMv is the unique tensor defined by hfivu" = and h/^vh'^^ = Sfj,^ — Tf^u^. If one 
has two sound cones, as in isotropic elasticity, it is the faster ray cone which lies 
outside. In Ex.3 C is given as a subset of vectors X'^ in a linear space T, which 
is the annihilator of the null space of s^", namely where this subset is given by 
s^vX^X'^ — 0, where s^^ is the inverse of s'"' on T. In the magnctohydrodynamic 
example the preceding statement corresponds to the fact that Alfven waves "travel 
along the direction of the magnetic field" . For Ex.4 the ray cone C is a 4th order 
cone of the same type as 6* , a fact already known by Ampere in the case of crystal 
optics and shown generally in [36]. For anisotropic elasticity the structure of the 
ray cone does not seem to be fully known, except for a general upper bound on its 
degree, namely 150 on groimds of general algebraic geometry (see [15], [37], [21]) 
and detailed studies for certain specific crystal symmetries - which give rise to a 
beautiful variety of acoustic phenomena [44]^. 



3 Initial Value Problem 

We now come to the issue of posing an initial value problem for hyperbolic equations 
of the form of Eq.(l). This requires two things: firstly a notion of "spacelike" initial 
value surface, secondly a notion of domain of dependance. Not surprisingly these 
notions can be formulated purely in terms of the characteristic polynomial. A hy- 
persurface S in R" will be called spacelike, if it has a conormal lying everywhere 
in r*{l;) for some If the equation (1) is nonlinear, every property concerning the 
characteristic polynomial has to refer to the data of some reference field /o, i.e. 
the value of /o on S and those of its derivatives up to order Z — 1. It is then the 
case that E is spacelike also for any sufficiently near-by data. The reason is that 
^' e r*{^) can be characterized by Ai(^, > 0, and the eigenvalues Aj, being zeros 
of a polynomial having real roots only, depend continuously on the coefficients of 
this polynomial [1] . A point x in R" is said to lie in the domain of dependence of E if 
each causal curve (i.e. each curve whose tangent vector X satisfies {X, ^ for all 
^' G r*{£_)) through x which is inextendible intersects E exactly once. The Cauchy 
problem for Eq.(l) is said to be well-posed if, for the above data, there is a unique 
solution in some domain of dependence of E and, secondly, if this solution depends 
in some appropriate sense continuously on the data. The question then is whether 
well-posedncss holds under the above conditions. The answer is affirmative when 
Eq.(l) is linear with constant coefficients and the lower-order terms are absent. Then 
the initial value problem can be solved "explicitly" by using a fundamental solution 
("Green function" in the physics literature) - wliic;li in turn c;an b(^ obtained c;.g. 
by the Fourier transform. By a refined version of a well-known argument in physics 
texts concerning the wave equation in Minkowski space (see e.g. [3]), one can show 
that the fundamental solution is supported in r{^_), which is the closure of the set of 
causal vectors just described. The set r{£^) is a closed, convex cone, dual to r*{^). If 
the outermost component of the cone C(^) dual to C*(^) is convex, its closure is the 
same as /^(C), otherwise its convex closure is the same as r{(^). If one is interested 
in finer details than just wellposedness, even the linear, constant-coefficient case 

® There are computer codes designed for algebraic elimination, which might be worth 
applying to this problem [23] 
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becomes very nontrivial. An example is the question of the existence of "lacunas" , 
i.e. regions in r'(^) where the fundamental solution vanishes. For isotropic elasticity 
mentioned in Ex.4, when C2 < Ci (which is the experimentally relevant case), the 
fundamental solution vanishes inside the inner shear cone determined by C2. (Note 
that "inner" and "outer" are interchanged under transition between normal and ray 
cone.) For anisotropic elasticity this issue, or the somewhat related question of the 
detailed time decay, already presents great difficulties (sec [15,39]). The existence of 
lacunas for general, linear hyperbolic systems with constant coefficients was studied 
in [3]. 

The problem now is that many field equations in physics give rise to variable co- 
efficients, to various forms of lower-order terms and-or nonlincaritics. But if one 
has a system of PDE's with hyperbolic characteristic polynomial (such systems are 
often called "weakly hyperbolic"), which in addition has a well posed initial value 
problem, a perturbation of the coefficients will in general destroy the latter property 
(see e.g. [31]). It is thus hard to get any further without additional assumptions. 
One such assumption is that of having a symmetric hyperbolic system. This is given 
by a system of the form of Eq.(l) with Z = 1. It is furthermore assumed that 

M/b=M^1b) (11) 

and that there exists so that 

MAsiO = ?m is positive definite. (12) 

The symmetric hyperbolic system has a characteristic polynomial which is hyper- 
bolic with respect to ^. To see this one simply observes that the equation 

det{M/j,{r]^ + X^^))=0 (13) 

characterizes eigenvalues of the quadratic form MAsiv) relative to the metric 
Mab (0 " ^'^d these eigenvalues have to be real. There is then, for quasilinear sym- 
metric hyperbolic systems, a rigorous existence statement along the lines informally 
outlined at the beginning of this section [29] . The uniqueness part uses the concept 
of "lens-shaped domains" (see e.g. [18]) which is essentially equivalent to that of 
domain of dependence above. 

Several field theories of physical importance naturally give rise to a symmetric hy- 
perbolic system. An example is afforded by the hydrodynamics of a perfect fluid 
both nonrelativistically and relativistically ^. The most prominent examples are 
perhaps the Maxwell equations in vacuo and the vacuum Bianchi identities in the 
Einstein theory. For the latter this was flrst observed in [17]. For completeness we 
outline a proof for the well-known Maxwell case following [45] . We have that 

V'F^, = 0, V[^i^,A] = (14) 

with being the covariant derivative with respect to g^i^, a Lorentz metric on Mf^. 
These are 8 equations for the 6 unknowns F^y. Next pick a timelike vector field 
with = —1 and define electric and magnetic fields by 

E^, = F^^u", B^^x = 3i^[^.WA], (15) 

so that 

F^, = -2E^u,^ - B^^xu^. (16) 

We assume for simplicity that u'' is covariant constant, otherwise the ensuing equa- 
tions contain zero'th order terms which are of no concern to us. The operator 

For an elegant treatment of the latter, see [16] 
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V^, = 2m[^V„] (17) 

contains derivatives only in directions orthogonal to u'^. Using Eq.'s (14) we find 
the evolution equations 

3V[^,Sa] = -uPVpB^.x, V^^'B^xp = 2uP\/pE,. (18) 

Taking now u^Vp of Eq.(15), wc rewrite the evolution equations in the form 

W^^^'^VaF^v' = 0. (19) 

Now take the positive definite metric 

w"^ = 2uV + g^"". (20) 

Consider now the positive definite metric a'^"^'' = 2w''^'^w''^^ on the space of 2-forms 
and use it to raise indices in W^J^ ^: One obtains quantities W^"^ " ^ satisfying 

Thus the Eq.'s (18) are symmetric hyperbolic with respect to u^. For the charac- 
teristic polynomial one finds 

P{k) ~ {u^k^figP^kpKf. (22) 



We now turn to 2nd order equations. Let us assume that the quantities M'^^ of 
Eq.(l) satisfy 

MTb = M''A- (23) 

This is necessarily the case when Eq.(l) comes from a variational principle, because 
then ^ 

where C = C{x, f, df). Of course, since only the quantities M^j^^^ enter the dif- 
ferential equation (1), one might as well have assumed the stronger conditions 
M^^ = = M^^. This is not usually done in continuum mechanics. The reason 

is that, while the PDE (and hence the characteristic polynomial) is unaffected by 
the above symmetrization, other physical quantities in the theory, like the stress, 
depend on the unsymmetrized object - and such objects typically enter, if not the 
equation, then the natural boundary conditions for the equation on the surface say 
of an elastic body (see [8]). A related reason is that the symmetrized object would 
hide other symmetries - present in some situations - which are more fundamen- 
tal such as invariances under isometrics. For example the object Cjt^(B\D\E)j with 
Cabde the elasticity tensor of Eq.(7), is not symmetric in {AB) and {DE). But 
the latter symmetry is important for understanding the solutions of the linearized 
equations of motion when the spacetime has Killing vectors. The work [10] also uses 
the unsymmetrized form of , the reason being that in this approach one is only 
interested in properties of M'^g which do not change when a total divergence is 
added to the Lagrangian, and the stronger symmetry, if present, would in general 
be destroyed by such an addition. Next it is assumed that there exists a pair X^, 
satisfying 

M^^^^^y is negative definite (25) 

and 

M'^j^{m'^r}^){m^r}^) > for all m^ry^ ^ with (X, 7]) = 0. (26) 
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The conditions (25,26) essentially state that the PDE is the sum of a "timelike part" 
and an "elliptic part" , the latter obeying the Legendre-Hadamard condition of the 
calculus of variations [22]. If the equation (1) has I = 2 and satisfies (23,25,26), 
the system is called regular hyperbolic with respect to ^. We now check that every 
regular hyperbolic system with respect to ^ is weakly hyperbolic with respect to ^. 
The characteristic condition reads 



rfet(M^^(r?^ + A^^)(r?. + A^.)) = (27) 

The covector r] in Eq.(27) can be decomposed as rj = \x'2) ^ ~^ ^ where I satisfies 
{X, I) = 0. Thus we can after redefining A assume that 77 in Eq.(27) has {X, rf) = 0. 
Defining Gab = -Mab{£) = -M'XUi^i,^ Vab = Mab{v), Qab = M^'XB^if.^u, 
consider the eigenvalue problem 

■Df = X£f, (28) 

in 



f = 



U 

where the quadratic forms V, £ are given by 



Vab 
Vab ^Qab 



and 



P_ [Vab \ 
V Gab)- 



Since £ is positive definite, all eigenvalues A are real. But Eq.(28) for / 7^ is 
equivalent to 

{-GabX' + 2QabX + Vab)v^ =0, (29) 

for which in turn is equivalent to Eq.(27). This proves our assertion that 

regular hyperbolic systems are weakly hyperbolic. (Note that every "timelike vec- 
tor" X in the sense of Eq.(26) is causal, i.e. (X,^) ^ for all ^ G r*(C), but not 
conversely.) We can now come back to Ex.5. The leading-order coefficients 
in Eq.(4) clearly belong to a regular hyperbolic system, when we choose the vector 
Xf^ uf^_ It then follows from the preceding result that the polynomial in Eq.(6) 
is indeed a hyperbolic polynomial. 

As with symmetric hyperbolic systems, it turns out that there is, for regular hy- 
perbolic systems, a local existence theorem [27] along the lines sketched at the 
beginning of this section. The appropriate domain of dependence theorem is proved 
in [10]. 

One may ask the question if it is possible to convert a regular hyperbolic system 
into an equivalent symmetric hyperbolic one by introducing first derivatives as ad- 
ditional dependent variables (at the price of course of having to solve constraints 
for the initial data). (This was the approach we originally followed for elasticity 
in [6], since we were unaware that there was already an existence theorem which 
applied, namely [27]). If the condition (12) is provisionally ignored, it turns out this 
is possible provided that is of the form of Eq.(4) for some pair u^, t^, i.e. 

certain cross-terms vanish But the positivity condition (12) will not always be 
satisfied. (Essentially this requires the "rank-one positivity" condition Eq.(26) to 
be replaced by the stronger rank-two positivity: M^^m^^^m^^ > for all m"*^ ^ 

® In [7] 1 claimed this to be possible even without these cross-terms vanishing. 1 now see 
1 have no proof of this assertion. 
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with X^m^ij, = 0.) In the case of isotropic elasticity it was shown in [6] that one 
can add to a term of the form Aj^^, which has the symmetries Aj^^ = ^[^gj, 
so that both the field equations and the requirement Eq.(23) remains unchanged, 
but at the same time condition (12) is vahd. However it is an algebraic fact that 
such a trick does not always work (see [38,42]). 

Finally let us mention the notion of strong hyperbolicity, which is intermediate 
between weak hyperbolicity and symmetric or regular hyperbolicity in the first or 
second order case respectively. This notion, which involves the tool of pseudodiffer- 
ential reduction [40,41], also gives wellposedness but has greater flexibility, see [33] 
for applications to the Einstein equations. It would be interesting to see if the chain 
"weakly hyperbolic - strongly hyperbolic - symmetric or regular hyperbolic" has an 
analogue for PDE's of order greater than 2. 
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